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The metal-insulator transition in correlated electron systems, where electron states 
transform from itinerant to localized, has been one of the central themes of condensed 
matter physics for more than half a century. The persistence of this question has 
been a consequence both of the intricacy of the fundamental issues and the growing 
recognition of the complexities that arise in real materials, even when strong repul- 
sive interactions play the primary role. The initial concept of Mott was based on 
the relative importance of kinetic hopping (measured by the bandwidth) and on-site 
repulsion of electrons. Real materials, however, have many additional degrees of free- 
dom that, as is recently attracting note, give rise to a rich variety of scenarios for a 
"Mott transition." Here we report results for the classic correlated insulator MnO 
which reproduce a simultaneous moment collapse, volume collapse, and metallization 
transition near the observed pressure, and identify the mechanism as collapse of the 
magnetic moment due to increase of crystal field splitting, rather than to variation in 
the bandwidth. 

We consider, as one of the simpler examples of the canonical Mott insulators,-'- the rocksalt 
structm'e transition metal monoxide (TMMO) manganese oxide with half-filled 3d shell. MnO 
is, most certainly, a multiorbital multielectron system with the accompanying complexities of the 
tenfold degeneracy, but the half-filled 3d states under ambient conditions lead to a spherical spin- 
only magnetic moment. Applying pressure to such a system leads to a number of changes, including 
insulator-metal transition, orbital repopulation, moment reduction, and volume collapse if a first- 
order transition results. These changes may occur simultaneously, or sequentially over a range of 
volumes.'^ Any of these may be accompanied by a structural phase transition, that is, a change in 
crystal symmetry, but an isostructural volume collapse may occur as well. The 3d bandwidth of 
such a Mott insulator is very susceptible to applied pressure, and has been thought to be one of 
the main controlling factors in the transition. 

While MnO's half- filled shell seems to give it a connection to well studied models, this aspect 
also makes it atypical of transition metal monoxides, as shown by Saitoh et al. who compiled^ ef- 
fective parameters for TMMOs from spectroscopic information. An effective intra-atomic Coulomb 
repulsion energy Ueff as defined by them is roughly twice as large as for the other 3d monoxides, 
and this has been used to suggest that MnO may be the most strongly correlated TMMO. The 
complexity that should be expected can be grasped by listing the relevant energy scales: 3d band- 
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width W, Coulomb repulsion U, intra-atomic d — d exchange energy (Hund's rule) J, crystal field 
splitting = — £t2g 1 and charge transfer energy Act = £t2g ~ [the difference in Mn 3d (we use 
t2g) and O 2p site energies]. All of these scales evolve as the volume decreases, altering the various 
microscopic processes and making the pressure-driven Mott transition a challenging phenomenon 
to describe. 

Early shock data^, and then Raman and optical studiesj^*^ had identified a transformation in 
MnO in the neighborhood of 90-105 GPa. Transport^^ magnetic," structural and spectroscopic,—*^ 
and reflectivity^ data all point to a first-order, insulator-metal Mott transition near 100 GPa with 
(reduced) volume {v = V/Vo) collapse t;=0.68^0.63, and moment collapse (from ~5/ub to 1^_b or 
jggg9iio^_ The structural data indicates a Bl— >B8 change just before the Mott transition, which thus 
occurs within the B8 (NiAs) phase rather than the Bl (NaCl) phase. Since the local environment 
of the Mn ion remains the same, this structural change is not expected to have much effect on the 
Mott transition in the disordered phase. 

Dynamical mean field theory (DMFT)iiii^ii^ as an approach for studying real materials has 
been showing impressive successes . ^^i^^i^^i^'^i^^'-^^'^^'^^' ^ The method that we have implemented and 
applied (see Methods section below) moves the treatment significantly beyond the methods used 
earlier for TMMOs, by including a full thermodynamic average of local dynamic processes result- 
ing from the strong interaction and all orbitals that can be relevant. Cohen, Mazin, and Isaak 
calculated the energy and magnetic moment using only local density approximation (LDA) based 
interactions.— In LDA MnO metalizes at (much too) low pressure; within the metallic phase they 
obtained a moment and volume collapse around 150 GPa. Fang and collaborators addressed this 
difficulty by using LDA only for the high pressure phase, and modeling the low pressure phase with 
the correlated LDA+U method.— With two different functionals, however, it is not possible to ob- 
tain the transition pressure. Four correlated electronic structure methods^^, applied throughout the 
volume range of interest, have probed the behavior of MnO under pressure; all obtained a high spin 
(HS, 5 = |) to low spin (LS, S = ^) moment collapse but their predictions differed considerably 
in other respects, demonstrating that the specific treatment of correlation effects is crucial. The 
prediction of the LDA-I-U method, which is regarded as the static, T=0 limit of the LDA+DMFT 
theory used here, is found to be affected by magnetic order,^ and predicts a zero temperature 
moment collapse in an insulator-insulator transition around 120 GPa (the pressure depends on 
the value of J), with little difference between the Bl and B8 structure results. Thermodynamic 
fluctuations have not been included in any previous study of MnO. 
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FIG. 1: This figure shows how the decrease in the local moment correlates with the orbital occupations, 

which reveal the spin state with relative volume.. The upper panel gives average instantaneous local moment 

Ms (black) and effective local moment Mef / (green) , and comparison to the Mn 3d orbital occupancies (lower 

panel) resolved into Cg (red) and (black) components. The solid lines represent the results obtained with 

the physical values U=6.9 eV, J=0.86 eV; the dashed lines using the enhanced value J=l eV and constant 

U/ J ratio illustrates how the moment collapse is suppressed to smaller volume if the spin-exchange coupling 

. . 4 

is increased. Closing of the t2g and Cg gaps is indicated by the black and red vertical lines respectively, 
confirming a connection between metallization and moment collapse. The analogous closing of the gaps for 




Energy (eV) 

FIG. 2: Ambient pressure x-ray photoemission spectroscopy (XPS) and bremsstrahlung isochromat spec- 
troscopy (BIS) data of van Elp et aZ. 30] on both sides of the energy gap for MnO. The upper curve, offset 
for clarity, is to be compared with the present DMFT resuh (bottom curve). While the separation of the 
main peaks is underestimated by 10%, the overall agreement in positions of structure is excellent. 

Magnetic Moment Collapse and Metallization 

Following most closely the approach developed and implemented by McMahan, Held, and 
Scalettar^i^SiSg f^^. 

pressure studies of elemental lanthanides, we have addressed the pressure-driven 
collapse of the correlated insulating state, using MnO as the prototype. Fig. [1] illustrates the evo- 
lution of the local magnetic moment and Mn 3d occupancies with volume. We use two different 
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measures of the local moment: (a) the mean instantaneous moment defined as an equal time 
correlation function Mg = \/ {nil) (b) effective local moment defined through the local spin 
susceptibility Mg// = ^/Txioc- These two moments have similar T-independent values in materials 
with Curie- Weiss behavior. Under compression, the local moment and Mn 2>d orbital occupancies 
retain their ambient pressure HS values (5 = |) down to about f=0.68. Further compression 
rapidly degrades the moment, which is accompanied by redistribution of electrons eg t2g within 
the Mn 3d shell. The local moments and orbital occupancies start to level off to the LS values 
around v = 0.57. The reduction of M^fj below Mg in the LS state indicates that the local moment 
screening (charge fluctuations) is enhanced in comparison to the HS state. 

Next we address the spectral properties, where the shortcomings of the LDA spectrum have been 
clear for decades. In Fig. [2] we compare the calculated total Mn 3d spectral function at ambient 
pressure with the photoemission data of van Elp et al."^^. Excellent agreement is obtained for the 
gap and for the peak positions. (We note that using the enhanced value of J=l eV gives significantly 
poorer agreement.) Having obtained a correct ambient pressure spectrum, we proceed in the study 
of the Mott transition by following the evolution with decreasing volume of the symmetry-resolved 
(t2g, Gg) spectral densities, presented in Fig. [3l The onset of the moment collapse around v = 0.68 
is signaled by, and associated with, closing of the gap in the t2g channel, while the eg gap is still 
visible at v = 0.63. This orbital selectivit y^-^ in metalization cannot be an exact property since 
both eg and t2g bands hybridize with the same O 2p bands throughout the Brillouin zone; however, 
the smallness of t2g-2p mixing allows the orbital selectivity to be remarkably pronounced. As the 
t2g gap closes, a quasiparticle peak appears at the chemical potential (E=0) as has been seen in 
simple models. Once in the LS state, the spectral functions bear strong resemblance to the parent 
LDA bands. In particular, the LDA (U = J = 0) t2g spectrum contains a sharp peak just at/below 
the chemical potential, so it is not certain how much of the peak arising at the transition is due to 
the many-body nature of the system. 

Mechanism of the Mott Transition in MnO 

We now address a fundamental point of this work, namely the connection between moment 
collapse and metal-insulator transition, by observation of the impact of pressure on the effective 
Hamiltonian. Since U and J do not change, the pressure enters the calculation only through the 
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FIG. 3: View of the evolution of the Mn 3d spectral densities under pressure (pressure increasing from 
bottom to top). The single-particle spectral functions are resolved into eg (left) and t2g (right) irreducible 
representations for varying relative volume. Note the spectral weight shift under pressure: de-occupation 
of Bg occurs as the increase in occupation of t2g proceeds (occupation = integrated weight over negative 
energies). For the lowest volume we show the uncorrelated (LDA) spectra for comparison (dotted lines). 
Apparently the main spectral features at high pressure originate from the uncorrelated band structure with 
some many-body renormalization. At even higher pressures the spectra remain qualitatively unchanged with 
some reduction of the weight of the high energy shoulders. 

quadratic (one-electron) part of the effective Hamiltonian. Reducing the role of pressure down to 
fundamentals one ends up with two effects: (i) broadening of the 3d bands and (ii) increase of the 
crystal-field splitting A^/. (We define A^j in terms of the site energies of the Cg and t2g Wannier 
functions; the eg-t2g band sphtting is substantially larger due to ligand field effects.) The evolution 
of the leading band structure quantities, which are the nearest-neighbor hopping amplitude tpda-, 
Ac/, and A^, are shown in the inset of Fig. [5l Since the 2>d bandwidth arises mainly through Mn 
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3d ~ O 2p hybridization {W oc t^^/Act) the increase of tpd hopping with pressure is to some extent 
compensated by the overall lowering of the p bands (increase in Act). 

So far we have demonstrated a connection between the moment collapse and metal-insulator 
transition (MIT), yet the chicken-and-egg question - which property drives? which property fol- 
lows? - is not yet answered. To this end we have performed an additional calculation at v = 0.8 
(well within the insulating HS state) without any intra-atomic exchange (J=0). In spite of the 
large U and same U/W ratio, a LS solution is obtained, which is metallic although strongly renor- 
malized. This result clearly shows that the MIT is driven by the collapse of the moment, which 
cannot withstand the increase of Ac/. The transition is characterized as evolving from five half- 
filled bands t2g + eg (HS) to three t2g bands with one hole per site (LS). The interaction energy 
cost of moving an electron from site to site is determined by Ueff = + d^~^ — 2d^. Using 
the atomic configurations corresponding to HS and LS states one arrives at an effective repulsion^^ 
C/^^ = U + AJ = 10.3 eV and U^f'j- = U — J = 5.9 eV respectively, indicating much stronger 
inhibition of the electron propagation in HS state. Moreover, the Mn 3d{eg) - O 2p hybridiza- 
tion provides an additional screening channel for the effective t2g — t2g interaction in the LS state. 
Indeed, a calculation performed in the LS state with O 2p states integrated out (keeping the 3d 
bandwidth unchanged) before solving the interacting problem leads to more pronounced high en- 
ergy shoulders as compared to the solution with 0-p states fully included; an indication of stronger 
local correlations. 

These results establish that the transition is controlled by competition between the crystal-field 
spitting Ac/ (favoring the LS state) and the intra-atomic exchange coupling J (favoring the HS 
state). Although both energy scales are important for the outcome of the calculations, only the 
former (Ac/) is sensitive to an applied pressure. The importance of the value of J was also found 
in LDA+U studies of the Mott transitiouj^Si Recently Werner and Millis-^ studied a two band 
model with competing intra-atomic exchange and crystal-field splitting. In the parameter range 
relevant for the present study they found three different phases realized in the following order with 
increasing crystal-field splitting: (i) spin-polarized Mott insulator, (ii) metal with large orbital and 
spin fluctuations and (iii) orbitally polarized insulator. Moreover they found an orbitally selective 
closing of the gap upon doping in the vicinity to the (i)-(ii) phase boundary. The correspondence 
between their spin-polarized insulator phase (i) and the HS state of MnO is evident. The transition 
region in MnO and phase (ii) are both characterized by metallization and strong orbital fluctuations 
as well as the orbitally selective gap behavior. Also the LS state of MnO and the phase (iii) of the 
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FIG. 4: Schematic energy diagrams of the spin states at both ambient pressure and at high pressure in the 
collapsed phase. Occupations of the Mn 3c? levels are pictured for both high spin (left) and low spin (right). 
In the HS state two spin-up electrons occupy eg orbitals at the cost of 2Ac/ in energy, but the spin-exchange 
energy gain is -lOJ (5x4/2 — 10 pairs of parallel spin electrons). In the LS state, the crystal field energy 
cost has become too great, and although the spin-exchange energy is less [-4J from 3x2/2 (up) + 1 (down) 
— 4 pairs] there is a net energy gain. The LDA energy difference is also a factor. 

model exhibit similarity, the orbital polarization. The insulating character of their phase (iii) is 
dictated by band-filling and does not contradict the above analogy. 

Like almost all previous studies using LDA-I-DMFT, we have included only the density-density 
terms of the Coulomb repulsion. Although they are not expected to influence a first-order volume 
collapse (see the next section) especially above 1000 K, it is gratifying to obtain some confirmation. 
Werner and Millis used the full Coulomb interaction in their study, and the similarity of the behavior 
of their model to what is found here for MnO provides some verification of the unimportance of 
the neglected terms. 
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FIG. 5: The equation of state that quantifies the volume collapse transition is presented in two representa- 
tions. In the left panel is the internal energy versus volume (dashed line represents a spline interpolation), 
and the right panel shows the resulting volume versus pressure curve (shading indicates estimated uncer- 
tainty), obtained as a derivative of the spline interpolation of E(V). The red bar on the V{P) curve lying 
at the theoretical transition pressure P^*'* — 120 ± 15 GPa determines the volume collapse v — 0.68 0.59. 
The dotted curve represents V{P) for the enhanced value of exchange J=l eV, showing the shift of the Mott 
transition to higher pressure with larger J. Width of the shaded red bar indicates the uncertainty of P^:'^ 
due to stochastic nature of QMC. The inset in the left panel shows the evolution of selected tight-binding 
parameters (units of eV) ; note specifically the factor of three increase in Ac/ . 



Equation of State and the Critical Pressure 



To compare to high pressure experiments, knowledge of phase stability is needed, which can 
be obtained from free energy vs volume (equation of state EOS). The theoretical justification for 
applying DMFT using the underlying LDA description relies on a well-defined thermodynamic 
grand canonical potential functional, for which specific realizations have been suggested;i^i21i2^ 
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Since it is difficult to extract the entropic term^ ■ in the free energy we restrict ourselves to evaluation 
of the internal energy; in any case the variation of the entropy term is very small on the energy scale 
of several eV involved in the changes of total energy. We use the internal energy scheme of McMahan 
et a/ . 15,27,28,29 i^j-^^t is similar to that of Savrasov and Kotliar^^ corresponding to the expression 
E(y,T) = El£)a{V) + [EDMFTiV,T) — EMpiV)], where Elda is the all-electron (unpolarized) 
LDA energy, Ed met is the internal energy corresponding to the self-consistent (dynamic) DMFT 
solution for the effective Hamiltonian and Emf is the static mean-field internal energy. The EOS 
curve is shown in Fig. [5j The main feature is the deviation from convexity in the transition 
region, which leads to a calculated volume collapse v^^ = 0.68 0.59 at P^^ = 120 GPa. The 
metallization and moment collapse obtained here are not far from the high pressure data,— with 
the transition volume (pressure) being somewhat smaller (larger) than the experimental values 
^exp _ Q gg _^ 0.63, Pc^^ = 100 GPa. The state just above the collapse is a HS insulator, with the 
t2g gap about to close. The collapsed state is LS, the Cg gap having just closed making it metallic 
in both subshells. 



Summary and Outlook 



These results demonstrate that the underlying LDA band structure, buttressed by on-site inter- 
actions ([/, J) treated within the dynamical DMFT ansatz, provide a realistic description of the 
Mott transition in MnO without input from experiment. This study finally allows a determination 
of the mechanism of the transition, which could not be uncovered by experiment alone: the mag- 
netic moment collapse, volume collapse, and metal-insulator transitions occur simultaneously, but 
it is the increasing crystal field splitting (encroachment of the O^^ ion on the internal structure of 
the Mn ion) and not the increasing bandwidth that tips the balance. 

The current results illustrate success of the LDA+DMFT approach in describing a pressure- 
driven Mott transition in a strongly correlated insulator, joining the growing number of successes 
of this approach in other strongly correlated real materials. The Kondo volume collapse tran- 
sition in Ce^i^ and other elemental lanthanides,— and the realistic modeling of parts of the 
complex phase diagraroi^ and multiplet effects^^ in Pu reflect the progress in correlated metals, 
with low temperature properties (heavy fermion characteristics) remaining an imposing challenge. 
Impressive progress has also been demonstrated in the description of structurally-driver>i^i21 and 
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doping-driven'^^ metal-insulator transitions in transition metal oxides. Excitation spectr of 
the charge-transfer compound NiO at ambient pressure, where O 2p states are tangled with the 
3d states, have shown excellent agreement with experiment. These results on MnO bring an ad- 
ditional class of materials into the list of strongly correlated systems whose behavior is becoming 
understood due to recent theoretical developments. 

Theoretical Approach and Numerical Methods 

Single Particle Hamiltonian and Interaction Term 

The LDA+DMFT computational scheme^^ in its present implementation, applied previously to 
NiO,— proceeds in two steps: (i) construction of an effective multi-band Hubbard Hamilto- 
nian H via Wannier transformation from a converged (unpolarized, metallic) LDA solution cor- 
rected for double-counting of the on-site interaction, and (ii) self-consistent solution of the DMFT 
equations^i^ using the quantum Monte-Carlo impurity solver.— 

k,CT,a,/3 k,(T,7,5 k,(T,a,7 k,(T,7,Q! 

+ Y^ ^O'P ""iao-^-f/So-' ) 

where dko-o- (Pk7(7) is the Fourier transform of the operator di^a (pi-ycr), which annihilates the d (p) 
electron with orbital and spin indices aa (70") in the ith unit cell, and nf^^ is the corresponding d 
occupation number operator. 

The single particle part of the Hamiltonian was obtained by a Wannier function projection 
method^, which amounts to a unitary transformation in the Hilbert space containing Mn 3d, O 
2p bands and the next lowest empty (polarization) conduction band. The site energy of the Mn 3d 
orbitals was corrected for double counting of the d — d interaction by subtracting from the LDA site 
energy a Hartree-like term giving = — {N — l)lJnL£)A, where = 10 is the total number 
of orbitals per Mn site, U is the average Coulomb repulsion and ulda is the average occupancy 
per d-orbital. Since the p — d band separation Act in MnO, which is to some extent influenced by 
the choice of the double-counting term, is rather large in the transition region, small variation of 
Act will not change the results. 

The Coulomb Interaction Matrix 
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The on-site Coulomb interaction U^"^' {U, J) within the Mn 3d shell, restricted to density-density 
terms only, was expressed as usual^i in terms of the direct (U) and exchange (J) interaction 
strengths related to the Slater integrals -Fq) -^2, -^4- The numerical values^ of U=6.9 eV and J=0.86 
eV were obtained by the constrained LDA method^. Since they exhibit only a small pressure 
dependence, these values were used for all volumes. We used L=100 imaginary time points in the 
Monte-Carlo simulation, in which the standard single-field-flip moves were augmented by special 
global moves that played a crucial role in ensuring ergodic sampling in the transition region. To 
obtain an indication of the robustness of our results we perform, in parallel with these ab initio 
interaction strengths, calculations with an enhanced (by 15%) value of J=l eV (and fixed U/J 
ratio). All the presented results were obtained at the temperature T=1160 K, in the rocksalt 
structure. 

Monte Carlo Procedure; Introduction of Global Moves 

The DMFT equations were solved numerically on a Matsubara contour (using asymptotic expan- 
sions for frequencies uJn > 500 eV), and the k-space integrals were performed by summation over 
3375 k-points in the first Brillouin zone. The chemical potential was adjusted in each DMFT itera- 
tion to guarantee the total electron count of 11 it 10^^. The impurity problem was solved using the 
Hirsch-Fye QMC algorithm^ modified for multiple orbitals. The on-site interaction was decoupled 
using a single binary Hubbard-Stratonovich auxiliary field Sap{l) for each pair of orbitals a(3 and 
each of L imaginary time slices (45 auxiliary fields for each time slice). 

The key innovation in this application to MnO in the transition regime was introduction of global 
Monte-Carlo moves in addition to the usual single-flips of the auxiliary flelds. These moves allow 
for fluctuations between HS- and LS-like conflgurations, which are otherwise practically unreach- 
able with the standard single-auxiliary- field-flip moves. The purpose of global moves is to mimic 
transferring electrons between orbitals. In general there is no straightforward relationship between 
a given configuration of auxiliary flelds, described by a binary L-vector S^p and the occupancy of 
orbitals. However, in the case of two atomic orbitals the probability distribution is peaked around 
auxiliary field configurations corresponding to a particular orbital occupancy, and flipping all fields 
corresponds to swapping occupancies of the two orbitals^^. This can be generalized to multiple 
orbitals as follows. To swap occupancies of orbitals a and /5 one has to: (1) fiip fields in Safi, 
(2) for all remaining flelds coupled to orbitals a or /3 swap the configurations S^a S^p. Since 
the decoupling is anti-symmetric with respect to the ordering of orbitals, auxiliary fields must be 
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flipped in step (2) whenever the order of orbitals changes between S-ya and 5^/3. 

Testing several types of the above moves we found that only simultaneous moves of two electrons 
between t2g and Sg orbitals of opposite spin (i.e. moves intuitively expected in LS^HS fluctuations) 
have appreciable acceptance. The acceptance rate of the global moves was found to be large 
only in the transition regime, which had been characterized by unusually slow convergence of 
the DMFT cycle. We checked for the possibility of multiple solutions, but found none at the 
temperature of these simulations. The numerical value of the total energy, limited by the stochastic 
error of the E^mft term, was converged to the accuracy of 0.06 eV in the transition regime and 
0.02 eV anywhere else. The spectral densities were calculated by the maximum entropy analytic 
continuation technique^ applied to the imaginary-time Green functions from 4 x 10^ — 6 x 10^ 
QMC-simulation sweeps collected into 2000-20000 bins. 
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